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INTRODUCTION 

The  stresses  acting  on  a  buried  structure  due  to  the  effects  of  a  weapon 
blast  depend  on  the  type,  size, and  orientation  of  the  weapon,  the  distance 
between  the  weapon  and  structure,  and  the  stress-propagation  behavior  of  the 
surrounding  soil.  For  a  given  weapon  loading,  the  response  of  the  structure 
depends  on  the  stresses  in  the  surrounding  soil,  while  the  stresses  in  the 
soil  depend  on  the  structural  response.  This  coupling  of  effects  is  termed 
soil-structure  interaction.  Depending  on  the  mass  and  stiffness  of  the 
structure  and  the  behavior  of  the  surrounding  soil,  the  stresses  acting  on  the 
structure  may  be  greater  than  or  less  than  the  free-field  stresses  'stresses 
that  would  exist  if  the  structure  vaj  re  not  present). 

Current  Air  Force  design  procedures  (References  1-3)  recognize  the 
existence  of  soil-structure  interaction  effects,  but  a  rational  method  of 
including  these  effects  in  the  analysis  procedure  has  not  been  implemented. 

For  a  complete  soil-structure  interaction  analysis,  the  soil  and  structure 
should  be  analyzed  as  a  system  with  appropriate  material  constitutive  models. 

Unlike  the  dynamic  loads  from  machines,  earthquakes  or  water  waves,  blast 
loads  may  result  in  stresses  of  large  magnitude  with  short  rise  times  and 
duration.  As  the  stress  wave  propagates  through  the  soil,  the  stresses  are 
attenuated  and  the  rise  time  of  the  stress  history  increases  (or  the  loading 
rate  decreases).  Since  the  propagation  and  attenuation  of  the  stress  wave 
depends  on  the  Inelastic  stress-strain  behavior  of  the  surrounding  soil ,  and 
the  response  of  soil  may  depend  on  the  rate  at  which  the  load  is  applied,  the 
wave  propagation  and  loading  rate  effects  are  closely  linked.  Thus,  for  an 
accurate  numerical  analysis  of  the  soil-structure  system,  the  constitutive 
characterization  of  the  soil  is  of  utmost  importance. 

The  objective  of  this  research  is  to  investigate  the  wave-propagation  and 
constitutive  behavior  of  geologic  media  subjected  to  blast-induced  stresses 
from  conventional  weapons.  This  study  is  to  consist  of  the  following  three 
tasks : 
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1.  Review  previous  experiment.il  investigations  of  wave  propagation 
tlirough  geologic  media.  This  study  should  address  the  effect  of  the  inelastic 
stress-strain  response  on  the  stress  propagation  behavior. 

2.  Review  existing  constitutive  models  that  may  describe  nonlinear, 
inelastic,  rate-dependent  behavior.  Of  significant  importance  in  the  choice 
and  implementation  of  a  model  is  the  ability  to  determine  numerical  values  for 
t  lie  model  parameters  from  laboratory  tests. 

3.  Provide  recommendat ions  regarding  the  implementation  of  constitutive 
models  with  Loading  rate  effects  into  numerical  codes  available  to  the  Air 
Force.  These  recommendations  should  include  applications  of  numerical 
techniques  to  problems  of  interest  to  Air  Force  research  and  analysis 
personnel.  Suggestions  with  respect  to  additional  laboratory  testing  programs 
and  additional  instrumentation  in  currently  scheduled  field  tests  will  be 

of  tered . 

The  resulting  "state-of-the-art"  report  should  give  Air  Force  personnel  a 
better  understanding  of  the  wave  propagation  of  blast-induced  stresses  and  how 
Loading  rate  effects  may  be  included  in  analysis  procedures.  These  analytical 
methods  may  be  useful  when  performing  survivabi Lity /vulnerability  analyses 
and  in  providing  geotechnical  guidance  tor  the  siting  of  buried  and  partially 
buried  structures.  In  addition,  a  numerical  analysis  prior  to  the 
ti eld- testing  of  structures  may  be  lieiptui  in  optimizing  the  test  with  respect 
to  ins tr umentation  and  weapon  placement. 
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The  determination  of  the  treo-ttold  stresses  acting  against  a  buried 
structure  is  essentially  a  problem  of  wave  propagation  through  the  surrounding 
soil.  The  sudden  release  of  energy  by  a  detonated  weapon  results  in  a  stress 
wave  tint  propagates  through  the  soil  in  a  spherical  wave  front.  The  stress 
wave  is  geometrically  or  spatially  attenuated  as  it  travels  from  the  source 
since  tlK  detonation  energy  is  being  distributed  over  an  increasing  volume  of 
material.  In  addition  to  the  spatial  attenuation,  the  stress  wave  is 
attenuated  by  the  dissipative  properties  of  the  soil.  This  attenuation  effect 
is  a  material  property  which  is  included  in  the  constitutive  model  and  should 
not  be  contused  with  the  spatial  attenuation.  The  spatial  attenuation  will  not 
be  considered  turtlier  as  it  is  included  automatically  in  the  analysis 
procedure. 

Most  of  the  early  wave  propagation  and  dynamic  stress-strain  testing  was 
conducted  to  determine  the  effects  from  nuclear  air-blast  loading,  but  the 
loadings  due  to  conventional  weapons  are  similar  in  many  respects.  As 
described  by  Wilson  and  Sibley,  (Reference  4)  the  classical  air-blast  wave 
from  a  nuclear  weapon  results  in  a  loading  history  at  a  given  point  in  the 
soil  that  rises  sliarply  and  is  followed  by  an  exponential  decay  in  the  positve 
phase  of  the  pressure.  This  positive  phase  is  followed  by  a  negative  phase 
rrequently  ignored  in  analysis  procedures.  Since  the  radial  extent  of  the 
Loaded  area  is  large  with  respect  to  the  thickness  of  the  soil  layer,  the  soil 
is  laterally  constrained  and  t i»e  resulting  deformation  is  essentially 
one-dimensional.  The  assumption  of  one-dimensional  loading  may  not  be 
appropriate  lor  tit  analysis  of  conventional  weapon  effects  against  buried 
structures  because  of  tlw  spherical  wave  front  that  is  generated.  Also,  the 
stress-tune  history  of  a  nuclear  weapon  has  a  slightly  longer  rise  time  and 
lasts  longer  tlian  loadings  usual  Ly  associated  with  conventional  weapons.  In 
spite  ot  these  differences,  many  similarities  exist  and  an  understanding  of 
one-dimensional  wave  propagation  is  instructive  and  may  be  an  appropriate 
assumpt  ion  lor  some  loadings  and  geometries. 


Tilt1  sliape  ol  tilt-  stress-strain  relationship  lias  been  recognized  as  the  most 
important  iactor  at  tec  ting  the  wave  propagation  response  of  soils  (Reference 
4-7).  As  indicated  in  Figure  1(a),  in  an  infinite  column  of  an  elastic, 
homogeneous,  material,  a  one-dimensional  stress  wave  propagates  with  a 
velocity  proportional  to  t  lie  elastic  modulus,  fig,  with  little  attenuation  or 
change  in  sliape.  In  one-dimensional  compression,  many  soils  exhibit  the 
tamiliar  "s-shaped"  response.  Figure  1(b).  Since  the  unloading  modulus,  Mu , 
is  larger  tlian  the  loading  modulus,  Mj ,  the  unloading  wave  propagates  more 
rapidly,  interacting  with  the  loading  wave  and  resulting  in  an  attenuation 
with  distance  trom  the  source.  As  indicated  in  Figure  1(b),  the  form  of  the 
stress  pulse  changes  sliape  as  it  propagates  through  the  soil,  resulting  in  an 
increase  in  tlie  rise  time,  a  decrease  in  tlie  peak  stress,  and  an  increase  in 
t  lie  total  duration  of  tlie  at  low  stress  levels,  the  relatively 

large  value  of  initial  modulus,  Mg ,  Figure  2,  results  in  a  stress  wave 
propagation  at  about  tlie  seismic  velocity,  while  at  higher  stress  levels,  the 
smaller  value  of  Mg,  results  in  a  propagation  velocity  less  than  the  seismic 
velocity.  This  explains  tlie  high  velocities  obtained  from  seismic  surveys  in 
which  relatively  low  stress  levels  ace  applied  and  indicates  tlie  importance  of 
testing  ac  appropriate  stress  levels. 

A  review  of  wave  propagation  can  be  considered  to  include  two  related 
areas  of.  study.  First,  experimental  studies  of  tlie  propagation  and  attenuation 
of  stress  waves  through  soils  will  be  reviewed.  In  these  tests  a 
o ne-d imensional  compression  wave  is  applied  to  one  end  of  a  long  soil  specimen 
and  t  lie  stress  and/or  velocity  recorded  as  the  wave  propagates  along  the 
s |iec  imen.  Since  tit*  wave  propagation  depends  on  the  stress-strain  behavior 
of  tit  soil  in  which  tit  wave  is  propagating,  tit  second  review  area  will 
consist  oi  efforts  to  measure  tit  dynamic  stress-strain  response  directly.  To 
measure  this  response  a  shock  wave  is  applied  to  a  short  specimen.  The  shock 
wave  will  travel  through  tit  specimen  and  be  reflected  several  times  during 
tit  duration  of  tit  pulse,  resulting  in  essentially  a  uniform  stress  state. 
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The  propagation  of  air  shock  waves  through  cylindrical  specimens  of  soli 
has  been  reported  by  several  investigators  (Reference  S-iO).  The  common  features 
in  these  tests  are  the  radiaL  confinement  of  the  specimen  and  the  application 
ot  a  stress  pulse  which  propagates  one-dimensional Ly  along  the  specimen.  The 
radial  confinement  is  provided  by  rigid  rings  or  spiral  reinforcing  in  a  rubber 
membrane  such  that  the  longitudinal  transmission  of  the  stress  wave  is  not 
significantly  impeded.  The  stress  and  displacements  are  measured  at  various 
locations  aLong  the  specimen  by  embedded  stress  and/or  strain  gauges.  While 
questions  may  arise  as  to  the  accuracy  of  stress-strain  measurements  taken  in 
this  manner,  the  overaLL  wave  propagation  and  attenuation  response  can  be 
effectively  observed  as  functions  of  soil  type  and  input  wave  form.  The 
"s-shaped"  stress  strain  response  has  been  recorded  in  shock  tests  on  sand 
specimens  confined  by  spiraL  windings  (Reference  7).  It  was  concluded  that  the 
formation  of  a  steady  shock  wave  in  sand  is  influenced  by  the  characteristics 
of  the  stress-strain  curve,  and  that  the  attenuation  dependson  the  material 
properties  and  boundary  conditions. 

Measurements  of  the  pore  air  pressure  in  samples  of  sand,  silt  and  pea 
gravel  (Reference  8)  indicate  that  pore  air  pressure  is  affected  by  the  peak 
overpressure,  shape  of  the  shock  wave,  permeability  of  the  material,  and  depth 
of  penetration.  Significant  pore  pressures  can  develop  in  coarse-grained 
uniform  soils,  but  are  of  little  consequence  in  well-graded  soils.  The 
attenuation  rate  is  independent  of  the  overpressure  and  depends  on  the  pore 
size,  depth  of  penetration  and  positive  shock  duration. 

Shock  wave  tests  on  Ottawa  sand  confined  in  segmental  rings  (Reference  9) 
indicate  materials  that  exhibit  strain-softening  display  an  increase  in 
the  rise  time  and  decrease  in  the  magnitude  of  the  stress  and  particle 
veLocity  pulse.  Strain  hardening  materials  tend  to  "shock-up”  or  show  a 
decrease  in  the  rise  time  as  the  wave  propagates  through  the  soil. 

Samples  of  Ottawa  sand  and  Kaolin  clay  (Reference  10)  have  been  subjected 
to  shock  waves  and  the  stress-strain  curves  for  the  sand  found  to  be  of  the 
tamilar  "s-shape,"  while  those  of  the  clay  were  essentially  linear  during  the 
loading  phase.  As  might  be  expected,  the  attenuation  rate  in  the  clay  was 
found  to  be  greater  than  that  ot  the  sand,  due  to  the  greater  hysteresis  loss 
•md  lower  wave  velocity. 
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Shock  wave  studies  have  been  performed  on  soLL  samples  without  rigid  Lateral 
confinement  In  which  confining  stresses  up  to  15  psi  were  applied  by  vacuum. 
Triaxlal  specimens  of  Ottawa  sand  (Reference  11)  and  clay  (Reference  12)  2.8 
inches  in  diameter  by  64  inches  long  were  subjected  to  air  shock  waves.  The 
results  indicate  that  the  stress  wave  in  sand  undergoes  the  least  amount  of 
change  in  shape  for  high  densities,  high  confining  pressures,  and  small 
amplitudes  of  shock  pressure  with  long  duration.  Because  these  conditions 
approximate  the  conditions  under  which  an  elastic  wave  would  exist,  not  much 
change  is  expected.  With  an  increase  in  the  peak  shock  and  decrease  in 
duration,  the  attenuation  was  found  to  Increase,  with  corresponding  increases 
in  peak  shock  and  duration. 

A  good  review  of  wave  propagation  studies  is  presented  by  Hampton  et  al. 
(Reference  13).  A  series  of  tests  on  laterally  confined  specimens  of  Kaolin 
is  described  in  which  rise  time  of  the  load  was  reported  as  "instantaneous”  or 
about  1  msec.  Measurements  from  embedded  stress  and  strain  gauges  indicated 
"softening"  behavior  or  stress-strain  curves  that  were  concave  downward 
towards  the  strain  axis.  This  is  a  departure  from  the  "hardening"  behavior  or 
"s-shaped"  response  typical  of  uniaxial  strain.  These  results  were  discussed 
with  respect  to  uniaxial  strain  tests  conducted  with  the  WES  device  (Reference 
14)  and  the  differences  were  attributed  to  a  difference  in  boundary 
conditions,  sample  size,  method  of  confinement,  and  the  difference  in  loading 
rate.  The  effects  of  submillisecond  rise  times  on  the  stress-strain  response 
will  be  subsequently  discussed. 

The  increase  in  strength  of  soils  due  to  the  rapid  application  of  load 
lias  been  recognized  since  the  1840s.  The  effect  of  nuclear  weapons  loading 
provided  the  Impetus  for  many  experimental  studies  during  the  early  60s. 
Although  much  of  the  early  data  consist  of  only  the  rate  effects  on  the  soil 
strength  and  Initial  modulus,  this  provides  an  understanding  of  the  general 
behavior  ot  soils  under  impact  loads.  Of  additional  Interest  is  the  apparatus 
used  to  apply  the  impulsive  load  to  the  specimen,  since  the  form  of  the 
applied  stress  pulse  is  a  function  of  the  loading  apparatus.  Stress  histories 
with  a  short  rise  time  and  duration  typical  of  conventional  weapons  may  be 
■  Ilf  limit  to  obtain  with  some  test  devices,  but  the  stress  strain  response 
must  be  measured  over  the  range  of  anticipated  Loading  rates. 


bample  geometry  is  an  important  consideration  when  determining 
constitutive  parameters  since  the  Knowledge  of  the  complete  state  of  stress  is 
required  it  toe  results  are  to  be  properly  incorporated  into  a  general  model. 
The  propagation  and  ret  lection  of  stress  waves  in  the  sample  also  depend  on 
the  geometry. 

The  importance  of  sampLe  configuration,  specifically,  lateral  constraint, 
has  been  discussed  by  Wilson  and  Sibley  (Reference  4).  It  was  suggested  that 
the  value  of  the  constrained  modulus  Mc ,  is  bounded  below  by  triaxiai 
compression  tests  due  to  sampling  disturbance  and  the  lack  of  lateral 
constraint.  The  value  of  Mc  is  bounded  above  by  the  modulus  calculated  from 
the  seismic  velocity.  This  value  is  high  due  to  the  effects  of 
preconsoildat  ion  and  the  low  applied  stress  levels.  The  actual  value  of  Mc 
should  be  between  tlie  value  computed  by  a  hydrostatic  test  and  a  uniaxial 
strain  test.  A  series  of  uniaxial  strain,  hydrostatic  compression,  and 
triaxiai  stiear  tests  on  samples  of  tuft  were  reported  by  (ihrogott  (Reference 
lb).  The  relationship  between  dynamic  and  static  strength  was  found  to  depend 
on  the  state  ot  stress  in  the  sample,  although  it  was  not  clear  what  these 
effects  are.  it  is  clear  ttiat  loading  rate  effects,  as  well  as  general 
stress-strain  behavior,  depend  on  stress  path.  A  review  of  impact  tests  under 
various  stress  paths  follows. 

Triaxiai  tests  on  dry  and  undrained  saturated  sands  were  conducted  by 
Whitman  and  Mealy  (Reference  lo)  with  time  to  failure  ranging  from  5  msec  to  b 
minutes  and  the  results  indicated  little  increase  in  strength  for  dry  and 
dense  saturated  sands.  Loose  saturated  sands  tend  to  develop  lower  values  of 
excess  pore  pressure  when  loaded  rapid  ly, resui t ing  in  a  strength  increase  of 
about  4u  percent  for  a  decrease  in  failure  time  from  b  seconds  to  U.U25 
sec  oiids. 

A  similar  study  on  the  rate  effects  on  a  saturated  remolded  clay  was 
reported  by  Richardson  and  Whitman,  (Reference  17)  and  it  was  suggested  that 
the  rate  effects  are  due  to  the  internal  migration  of  pore  water  causing 
larger  excess  pore  pressures  at  low  rates. 


An  ex  to  ns  i  vo  so 1  i  oh  ol  tests  on  (loose  Lake  clay  (Keteroiiee  18)  at  various 
water  contents  was  conducted  with  a  pneumatic  device  producing  failure  times 
ranging  t  roin  2  msec  to  1  hour.  An  increase  in  dynamic  strength  was  observed 
tor  alt  water  contents,  but  a  smaller  strength  increase  was  found  for  samples 
with  low  water  contents.  The  increase  in  secant  modulus  at  1  percent  strain 
was  round  to  lie  more  dependent  on  strain  rate  and  water  content  than  was  the 
strength,  increases  in  dynamic  strength  and  modulus  have  also  been  reported 
by  other  researchers  (Reference  19-21). 

Dynamic  direct  shear  tests  (Reference  22  and  23)  using  compressed  air  to 
create  loading  rise  times  ranging  from  1  to  5  msec  have  been  conducted  on 
specimens  of  sand,  clay,  and  sand-clay  mixtures.  The  dynamic  impact  strength 
envelopes  were  compared  with  those  of  "rapid  static"  tests  with  time  to 
failure  of  30  to  30  seconds.  The  results  indicate  little  increase  in  strength 
tor  tit;  sand  samples  in  both  dry  and  saturated  conditions,  while  the  clay  and 
clay-sand  mixtures  exhibit  an  increase  in  apparent  cohesion.  The  increase  in 
cohesion  appears  to  be  independent  of  the  moisture  content  and  grain  size  for 
specimens  with  a  degree  of  saturation,  S,  greater  than  85  percent. 

The  split  Hopkinson  pressure  bar  (Reference  24)  or  Koisky  (Reference  25) 
apparatus  was  developed  for  tic  testing  of  metals  and  has  been  used  to 
investigate  tic  loading  rate  effects  on  the  behavior  of  geologic  materials 
(Reference  26-29).  Although  the  Hopkinson  bar  is  convenient  for  dynamic 
testing  and  tic  impact  velocities  can  be  accurately  controlled,  the  rate  of 
strain  is  not  constant  and  cannot  be  determined  directly.  This  apparatus  and 
the  analysis  of  tic  test  data  have  been  well  documented  (Reference  30-33)  and 
will  not  be  discussed  further. 

Perhaps  tic  best  data  on  loading  rate  effects  have  been  obtained  from 
tests  conducted  on  thin  specimens  in  uniaxial  strain  devices  (Reference 
5,14,34,35).  Since  the  sample  length  is  small,  the  stress  wave  may  propagate 
across  tic  sample  and  be  reflected  from  the  ends  many  times  during  the  duration 
of  the  load.  Thus,  a  relatively  uniform  stress  distribution  throughout  the 
sample  can  be  assumed.  This  assumption,  known  as  the  "multipie-reflection 
theory,"  is  probably  appropriate,  provided  the  rise  time  is  greater  than  0.1 
msec  tor  sand  and  1.0  msec  for  clay,  (Reference  36)  but  reflections  and 
rarefactions  may  atlect  tic  actual  state  of  stress  in  the  sample  (Reference 
3/).  In  addition,  the  complete  state  of  stress  is  unknown  due  to  the  inability 
to  measure  tic  radial  stresses. 
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Schindler  (Reference  14)  describes  the  development  of  a  uniaxial  device 
by  WbS  in  which  a  cold-gas  expansion  is  used  to  apply  the  impact  loading. 

This  device  was  moditied  by  Jackson  et  al.  (Reference  35)  to  duplicate  the 
eftects  of  a  conventional  weapon  by  using  a  loading  initiated  by  a  hot  gas 
expansion  (explosion)  resulting  in  rise  times  less  than  1  msec.  The 
stress-strain  responses  of  three  sands  liave  been  measured  with  this  device 
under  transient  pressures  ranging  from  10  mPa  to  80  mPa.  It  is  reported  that 
the  response  changes  from  "stiffening"  behavior  to  "softening"  behavior  when 
the  rise  time  is  less  than  1  msec.  A  numerical  analysis  was  performed  to  show 
tint  these  eftects  were  due  to  rate  eftects  rather  than  a  change  in  the  state 
of  stress  due  to  reflections.  This  distinct  change  in  behavior  with 
subrail lisecond  rise  times  is  similar  to  that  reported  by  Hampton  (Reference 
13).  It  lias  been  suggested  (Reference  38)  that  this  softening  is  due  to  a 
decrease  in  strain  rate  after  the  initial  rapid  loading.  This  change  in 
behavior  is  significant  and  should  be  investigated  more  closely  under 
different  states  of  stress.  Further  testing  with  rise  times  less  than  l  msec 
should  be  conducted  to  assure  that  this  measured  response  is,  in  fact, 
representative  of  the  material  and  not  a  function  of  the  testing  device. 

Table  1  summarizes  the  rate  effects  on  strength  and  initial  modulus. 
Although  a  complete  description  of  the  stress-strain  response  is  needed  for  a 
wave-propagation  analysis,  the  general  trends  are  informative.  The 
stress-strain  behavior  of  the  soil  is  highly  dependent  on  the  type  of  soil, 
lateral  confinement,  state  of  stress  or  stress  path,  stress  level,  and  the 
rate  at  which  the  load  is  applied.  Laboratory  tests  should  be  conducted  under 
conditions  approximating  those  in  situ.  Some  existing  mathematical  models  that 
represent  the  constitutive  response  of  geologic  materials  will  now  be 
discussed.  Of  particular  interest  are  models  which  incorporate  the  effects  of 


TABLE  1.  .SUMMAKy  Oh'  EFFECTS  ON  STRENGTH  ANO  INITIAL  MODULUS  DUE  TO 
INCREASE  IN  LOADING  RATE* 


SOIL  TYPE 

STATE 

STRENGTH 

INCREASE 

INITIAL  MODULUS 

INCREASE 

COMMENTS 

Ail  S 

less  than 
increase  in 
modulus 

.  .  _  J 

great 

lower  excess 
pore  pressures 

function  of  OCR 

CLAY 

Low  S 

moderate 

function  of 
moisture  content 

S  "  U5Z 

great 

independent  of 
moisture  content 
and  grain  size 

SAND 

Dry 

rggn 

si ight 

slight 

dependent  on 
density 

6, 

GRaVKL 

Sa  tu  ra  ted 

loose 

great 

lower  excess 
pore  pressures 

de  nse 

slight 

slight 

*NOT E:  General  trends  are  indicated  only.  The  effect  of  loading  rate  is  a 

function  of  many  variables,  such  as  the  applied  stress  level,  the  stress 
state,  tin;  loading  History,  the  grain  size  distribution,  density,  water 
content,  porosity,  and  overconsoi  idat  ion  ratio,  amoung  others.  The  properties 
of  a  given  soil  should  be  determined  in  tlte  field  or  laboratory  under  ia  situ 
conditions,  and  shouid  be  determined  under  various  states  of  stress  and  over 
the  range  of  expected  loading  rates. 
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Sfc.cn.uN  in 


CONS  l'l  l  UTlVK  KtLATlONS 

Numerical  wave  propagation  (Reference  3V-41)  and  soil-structure 
interaction  modelling  tor  blast  loadings  (Reference  42-43)  nave  been  conducted 
witn  elastic  material  representations  in  which  the  stress  is  related  directly 
to  toe  strain.  Due  to  t lie  inelastic  response  and  t lie  stress-path  dependence  of 
most  soiLs,  an  incremental  plasticity  formulation  such  as  that  conducted  by 
Goudreau  et  al.  (Reference  44)  tor  a  combined  air  blast  and  ground  shock 
loading  is  preferred. 

Incremental  plasticity  type  models,  in  terms  of  stress  invariants,  can 
represent  material  behavior  over  multiple-stress  paths  and  will  assure 
uniqueness  of  tlie  stress-strain  relation.  Space  does  not  permit  a  complete 
discussion  of  plasticity,  but,  in  most  formulations  used  to  represent  geologic 
materials,  the  strain  is  assumed  to  consist  of  an  elastic  part,  ee,  and  a 
plastic  part,  eP,  as  siiown  in  Figure  2.  To  include  the  frictional  response 
common  to  most  geologic  media,  a  failure  envelope  such  as  tliat  shown  in  Figure 
3  is  used.  When  tlie  stress  state  (plotted  in  an  invariant  space  as  shown  in 
Figure  3)  is  below  the  envelope,  only  elastic  strains  can  occur.  Since  most 
geologic  materials  experience  plastic  strains  at  stresses  below  the  linear 
envelope,  or  due  to  stress  paths,  as  shown  In  Figure  3,  an  additional  yield 
criterion  or  cap,  is  defined.  As  the  stress  state  moves  from  Point  l  to  Point 
2  the  cap  also  moves  as  a  result  of  strain  hardening.  This  results  in  tlie 
change  in  tne  elastic  limit  of  the  material  response.  While  many  plasticity 
tormuiat  ions  exist,  (Reference  4b  and  4b)  formulations  that  can  represent  the 
loading  rate  effects  are  of  interest  In  blast  loadings.  A  trade-off  must 
often  he  made  between  sophistication  of  tlie  model  and  tlie  ability  to  determine 
tlie  material  parameters  from  laboratory  or  field  tests. 

Although  many  of  tlie  elas to-v lscop las 1 1c  models  are  formulated  to  model 
creep  behavior,  t'ie  formulation  sliould  be  general  enough  to  include  the  change 
in  behavior  due  to  rapidly  applied  loadings.  Variations  of  tne  Sandler  Cap 
model  (Reference  4d )  have  been  proposed  (Reference  3b,  49,50)  in  which  tlie  flow 
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rule  is  written  in  terras  of  a  viscous  parameter.  Under  slow  rates  of  loading, 
t  l»e  cap  moves  witli  tie  stress  state  as  in  a  conventional  plasticity  model,  and 
tie  hardening  or  stittening  behavior  expected  for  a  uniaxial  stress  state  may 
occur.  Under  rapid  rates  of  loading,  tie  cap  lags  the  stress  state  indicated 
by  Point  3  producing  the  "apparent  softening"  (Reference  49)  behavior  measured 
by  Jackson  (Reference  35).  Additional  strains  accumulate  in  the  visco-plastic 
domain  until  tie  cap  "catches  up"  to  the  stress  state.  A  model  of  this  type 
has  been  formulated  by  Rubin  and  Sandler  (Reference  49)  for  high-pressure  and 
high-temperature  equations  of  state,  and  modular  subroutines  are  provided  for 
use  in  explicit  integration  codes.  A  model  formulated  in  terms  of  the  third 
stress  invariant  has  been  used  by  Schreyer  (Reference  51).  Other  models 
(Reference  38,  49-51)  have  all  been  used  to  predict  the  laboratory  tests  of 
Jackson  (Reference  35)  and  appear  to  represent  the  rate  effects  adequately.  It 
should  be  kept  in  mind  that  due  to  the  lack  of  experimental  data,  only  the 
behavicr  under  the  uniaxial  state  of  stress  has  been  predicted  for  rise  times 
less  than  1  msec.  A  true  test  of  these  models  awaits  laboratory  data  under 
different  stress  paths. 


stcrioN  iv 


RECOMMKN  DAT  IONS 

Existing  numerical  techniques  are  available  to  perform  dynamic 
soil-structure  interaction  analyses  for  weapons  effects  on  buried  structures. 
For  many  soil  conditions  such  as  clays  and  partially  saturated  sands,  a 
rate-dependent  constitutive  raodeL  may  be  required.  Since  an  explicit  time 
integration  scheme  is  usually  recommended  for  blast  loadings,  (Reference 
31,52)  a  large  deformation  code  such  as  DYNA3D  (Reference  53)  may  yield  good 
results. 

To  indicate  how  such  a  code  may  be  used,  the  full-scale  buried  HYPAR 
structure  tests  to  be  conducted  at  Tyndall  AFB  should  be  analyzed.  Soil 
samples  from  the  test  site  should  be  tested  under  quasi-static  loading 
conditions  and  the  analysis  performed  with  constitutive  models  currently 
implemented  in  DYNA3D.  The  ground  motions  and  interface  pressures  will  be 
compared  with  those  measured  in  the  Held  tests  to  illustrate  the  importance 
of  the  soil  model.  While  a  better  solution  is  likely  to  be  obtained  if  a 
rate-dependent  model  is  used,  the  analysis  will  demonstrate  the  laboratory 
determination  of  the  rate  independent  constitutive  parameters  and  the  analysis 
with  soil-structure  interaction  effects.  With  rate-dependent  properties  from 
future  tests  at  high  loading  rates,  a  better  solution  may  be  obtained. 

Future  Air  Force  research  should  be  directed  towards  development  and 
implementation  of  rate-dependent  constitutive  models.  The  severe  shortage  of 
quality  laboratory  test  data  under  high  loading  rates  hinders  this-* 
development,  therefore  both  tasks  should  be  considered  concurrently.  Stress 
states  otlier  than  uniaxial  strain  are  important  for  the  localized  stress 
conditions  existing  near  conventional  weapons.  Data  on  soils  tor  which 
tul  1-scale  structural  tests  have  been  conducted  or  are  planned  should  be 
considered  to  verity  the  models.  Numerical  analyses  of  field  test  structures 
prior  to  testing  can  permit  parametric  studies  of  test  variables  such  as 
weapon  location  and  orientation,  as  well  as  assist  in  the  determination  of 
instrumentation  location.  A  rational  approach  to  the  design  of  field  tests 
will  permit  more  effective  use  of  test  structures  and  weapons  to  maximize  the 


d.iiu  out  a  mi'll .  testing  can  then  follow  to  determine  it  the  structure  behaves 
as  predicted.  Numerical  analyses  can  also  be  useful  in  isolating  the  effects 
of  stress  wave  reflections  in  laboratory  and  centrifuge  tests  at  high  loading 
rates . 

Additional  areas  of  research  in  soil  structure  interaction  research  include 
tlie  development  of  appropriate  interface  models  to  permit  relative  motions 
between  the  soil  and  structure.  These  relative  motions  may  be  particularly 
important  tor  surface  and  partiaiiy  buried  structures.  Although  interface 
models  exist,  these  models  cannot  include  the  effects  of  soil  arching.  Soil 
arching  is  the  tendency  tor  the  soil  to  transfer  stress  from  the  more  flexible 
portions  of  tlte  structure  to  tlie  portions  of  greater  stiffness.  Tests 
simulating  nuclear  weapons  have  shown  that  current  analysis  methods 
underpredict  the  strength  of  the  structure  due  to  these  effects  (Reference  3). 
Arching  ettects  may  be  enhanced  witn  the  rapid  application  of  load  associated 
witu  conventional  weapons.  Special  interface  elements  can  be  developed  to 
include  tliese  ettects,  witn  the  constitutive  parameters  determined  from  field 
or  laboratory  tests. 

The  analysis  of  buried  structures  subjected  to  conventional  weapon  loadings 
is  complex,  but  existing  methods  are  available  to  obtain  reasonable  solutions. 
Witu  tlie  nnplemontat  ion  of  rate-dependent  material  models,  the  quality  of  the 
solutions  sliouid  improve.  Numerical  solution  procedures  interpreted  in 
conjuction  with  tlie  results  of  centrifuge  and  field  tests  will  result  in  an 
improved  understanding  of  tlie  soil  structure  interaction  and  wave  propagation 
effects  in  these  problems. 
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